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ABSTRACT 

The problem of the calculation of ambient densities from SM-C 
NACE data is considered. A brief description is given of the data 
measurement method, followed by a description of both the theoret- 
ical and experimental data curves. Geometry, electrometer dis- 
tortions, and noise effects are then studied in terms of their effects 
on the ideal data form. From these considerations two data reduc- 
tion methods are evolved. The flrstis an iterative integration tech- 
nique that exploits the symmetry of the experimental data about the 
minimum angle of attack. For the analysis of geometry effects, a 
second method using interval averaging was developed and studied. 
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SYMBOLS 


Definition of terms used in this paper: 


Symbol 


A 

n 



e 


i v 0 1 

I Vi l 


s 


s 


o 


F(s) 

G(s) 

<P 

T 

oj 

At 


Definition 


Experiment orifice normal 
Satellite velocity 


Angle of attack 


arcos 


A 

n 


• Vsat 


L IVsatl J 

Most probable speed of a species of molecules with a Maxwellian 
distribution outside of the orifice 


Most probable speed of a species of molecules with a Maxwellian 
distribution inside the orifice 


The speed ratio 


|Vsat|cos0 

|V C I 


J IVsatl 

Maximum speed ratio = — — 

IV 0 I . 

The thermal transpiration function, the idealized ratio of source to 
ambient density 

Correction to F(s) that accounts for the experiment chamber’s 
geometry 

Geometry corrected thermal transpiration function = F(s)r? So (0) 
Roll angle of orifice normal 
Satellite spin period 

Angular spin frequency of satellite = 2n/T 
Data spacing in seconds 
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AN INVESTIGATION OF GEOMETRY AND NOISE CORRECTIONS 
TO SAN MARCO-C NEUTRAL ATMOSPHERIC 
COMPOSITION EXPERIMENT DATA 


I. INTRODUCTION 

The goal of the San Marco-C Neutral Atmospheric Composition Experiment 
(SM-C NACE) is the measurement of the primary constituents of the neutral at- 
mosphere. These components are:' He, O, N 2 , 0 2 > and Ar, with provision 
for the measurement of contaminant gases (CO, H 2 0, C0 2 ) having been included. 


The measurement process employs a double focusing, magnetic sector type 
spectrometer, with which the incoming gas beam, singly ionized by electronbom- 
bardment, is separated into its component gases by their charge-to-mass ratios. 
These source ion currents are then converted, by means of electrometers, to 
voltages for telemetry. These electrometers are of a linear range switching 
type. Knowing the instrument transfer function allows one to convert the telem- 
etry voltages back to source currents. These source currents combined with a 
knowledge of current-pressure calibrations then determine the mass spectrom- 
eter ion-source densities which are related to the ambient neutral particle den- 
sities. All of the electronic transfer functions and gas calibrations are essen- 
tially linear in nature. 

In terms of geometry, the orifice normal of the enclosed ion source is in the 
satellite’s spin equatorial plane to a good approximation (±2.0°). 

The San Marco-C satellite had a nominal orbit with initial apogee of 713 km 
and perigee of 213 km with a nominal spin period, T, near 8 seconds. 

The ideal source density form, F(s), is shown in Figure 1 plotted vs. angle 
of attack. The actual data contained noise of varying severity. Figure 2 shows 
the spin modulated telemetry voltages for N 2 vs. time. This example was taken 
from strip charts of the data. Figure 3 shows the geometric relation between 
angle of attack and roll angle. 

The theoretical relation between source and ambient densities can be deduced 
from elementary kinetic theory. If one assumes: 

(1) Free molecular flow (mean free path» chamber dimensions and tra- 
jectories are straight lines) 

(2) Instantaneous accommodation at the orifice plane in the instrument 
chamber (no geometry effects) 
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(3) Dynamic equilibrium in the gas flows 

(4) No surface reactions and/or outgassing, and, 

(5) Both source and ambient particles may be described by Maxwell- 
Boltzmann statistics; one arrives at the equation 

NjO) y.(j) = N 0 <j) V o 0) F(sO)) (1-a) 

where 

= jth component of source density 
N 0 *J) = jth component of ambient density 
VjG) = most probable speed of jth source component 
V 0 ^ = most probable speed of jth ambient component 
s<j) = V sat ' n/|V 0 ( -i)| , the speed 

A A 

ratio of the jth component where n is the orifice normal and (V sat • n ) = |V sat 
I cos 6 , where 6 is the angle of attack, and finally; 

F(s) = exp(-s 2 ) + s \/lr (1 + erf s) ( 1 -b) 

The physical content of this equation is that since the satellite travels through 
the surrounding atmosphere at velocities greater than mean ambient thermal ones, 
the chamber orifice at small angles of attack ( 10 1 < 90°) overtakes almost all 
of the gas. In contrast, at large angles (|0| > 90°) the ambient atmosphere has 
only a very small proportion of it's population that can overtake the satellite and 
hence enter the orifice. 

There are two limiting cases of interest for which F(s) assumes a rather sim- 
ple form. First, for s large and positive (i.e. s>l), F(s) is well represented 
by the linear approximation: 

F(s) ~ s (2) 

Secondly, for s large and negative (i.e. s < -1): 

F(s) ~ exp(-s 2 ) (3) 


2 



■> •> 

and thus is exponential. Thus, for a given s 0 = |V sat | /V 0> s = s 0 cos 6 and one 
may compute ranges of 6 for which these approximations hold. Both of these cases 
are of interest in the analysis of atmospheric data, however, it is the linear s 
region that is of primary importance in density calculations. In this linear s re- 
gion the ambient-source density equation (1-a) simplifies to: 

Nj Vj = 2 % /¥N o |V sat |cos0 (4) 

which is independent of ambient temperature . 


II. CONSIDERATIONS OF THE GEOMETRY CORRECTION TO THE 
THEORETICAL F(s) 

If one includes geometry effects, the spin modulation of the source density 
is no longer adequately described by F(s) for angles of attack, d , < 35 . For 
then, at near full ram, the directed motion of the particles becomes important. 
This leads to a modification of the spin modulation that may be described as: 

G(s) = n So (0) F(s) (5) 

where the subscript s 0 indicates that r? depends both on the angle of attack, 6 , 
as well as the speed ratio s 0 at 6 = 0, where t? Sq ( 0) > 1. The theoretical pre- 
diction of this modification is shown in Figure 4. The dotted line visible for small 
6 shows the uncorrected function. It can be seen from this that for s 0 = 10.0521 
a value characteristic of N 2 at 1150° K and I V sa ^ I = 8.3km/sec, the maximum 
effect is ~ 40% for 6 - 0 (full ram). The Vs 0 ( d ) used is that Predicted by 
Dr. J. C. Pearl of Goddard Space Flight Center. As tj Sq (d) is not an analytical 
function of 0 it has been digitized from graphs and then spline interpolated for 
use at all angles of attack. Note that r? So (0) is a monotonic decreasing function 
of increasing |0|. 

With the inclusion of this correction the principle effect remaining is that of 
surface-gas interactions. These effects as predicted by Pearl results in a con- 
stant linear scaling factor to F(s) for |0| > 35° and tj So = 77 . ( e > 0 ) for |0| 

< 35° , where e is the recombination coefficient. For the relatively inert N 2 
and the inert Ar, and He which are considered here, r? = 7? $o (0) > 1 and sur- 
face effects are neglected. 

Thus given the geometry corrected transpiration function, G(s), one needs 
to consider the effects of the finite bandwidth of the electrometers through which 
the ion-source currents are converted to telemetry voltages. These effects will 
determine the shape of the theoretically observable G(s). For SM-C NACE these 
bandwidths ranged from 10.0 to 1.3 Hz. The satellite spin frequency was a 
nominal 0.125 Hz. 
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IE. ELECTROMETER BANDWIDTH CONSIDERATIONS 

To consider the electrometer effects on G(s), i.e. on the source data, one 
needs both a model for the electrometer and a suitable representation of G(s). 
G(s) was represented as a Fourier series: 


M 

G(s) %? 


a N (s)cos(N0) + b N (s)sin(N0) + a 0 


( 6 ) 


The coefficients of this series, a N (s) and b N (s) were generated by first obtaining 
an analytical representation of G(s) = F(s) r) So (0) by spline interpolating r? So (0). 
Then, the coefficients were computed from the standard integrals using a 32 point 
Gaussian Quadrature routine. The value of M at which the series was truncated 
was determined by that M which minimized: 

2 [G(s) - G Fourier (s)] 2 (7) 

where the sum is over all points in the spin period. 

The Fourier series constructed in this manner was found to agree with the 
ideal G(s) to within 1% in the range - 90 < 0 <90 degrees. 

With this G(s) representation a model electrometer was constructed from an 
RC circuit approximation in the form of a low pass filter. The attenuation of the 
Nth harmonic was then given by : 


a N 

7n = l./Vl. + (N/Tf 2 ) 2 


( 8 ) 


and the phase shift by: 


0 N = arctan (N/Tf 2 ) (9) 

where f 2 is the electrometer bandwidth and T is the satellite spin period. Then 


G'(s) 


M 

2 a'w(s) cos (N0 + 0 N ) + b' N (s) sin (N0 + 0 N ) + a 0 
N = 1 N 


( 10 ) 


\ 
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for the post-electrometer transpiration function. The relation between the 
primed and unprimed coefficients is given by: 


/ a " N \ = /cos0 N sin0 N \ /a N \ 
\ b N/ 7N \-sin0 N cos 0 N / V b N/ 


(ID 


which is easily inverted as the matrix is orthogonal, to give: 

/ a N \ = /cos0 N -sin tf> N \ /a' N \ 

V b N/ Tn \sin^ N COS 0 N / \b’ N / 


( 12 ) 


This allows a transformation of the data to the pre-electrometer side if so 
desired. 


The theoretically observable G’(s) that results from the transfer of this Fou- 
rier representation through the model electrometer has the following 
characteristics : 

(1) In the region where s > 1, G(s) was found to be shifted, to a good ap- 
proximation, by the phase shift of the fundamental harmonic with a very 
small attenuation. In the worst case, for a bandwidth of 1. 3 Hz, and 

|0| > 35° this attenuation was ~ 0.5% and the phase shift 5.49°. 

(2) For |0| <; 35° the distortion of the geometry correction was more se- 
vere as the information is contained in the higher harmonics. Although 
the attenuation is at most 5% of the total curve amplitude for 1 . 3 Hz, the 
correction factor was attenuated by about 30%. For large bandwidths 
(>8 Hz) the distortion is ~ 0. 5%. 

Figure 5 shows the pre- and post-electrometer G(s) for s 0 = 10.052 and a 
bandwidth of 1.3 Hz. The essential linearity of the electrometer transformation 
is shown in Figure 6. If one regards the 0.5% attenuation as negligible, the elec- 
trometer transformation in the linear s region for |0| > 35° is merely shifting 

the roll angle axis by a constant. The transformation can thus be regarded as 
linear. 


IV. NOISE SMOOTHING: INTEGRATION OF G(s) TO DETERMINE DENSITIES 
AND ROLL PHASE ANGLES 

The electrometer effect has been shown to be well approximated by a linear 
transformation involving a change of phase in roll angle, for all bandwidths, with 
the restraints that s > land |0| > 35°. The condition s > 1 implies |0| 

< 75° for values of s 0 as small as 4, which is characteristic of He at 600° K. 
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Thus, for 35° < |0| < 75°: 


G(s) = 2y/lt s 0 cos0. (13) 

From (2) and (4) the spin modulated source density is given by: 

/N 0 |V sat l\ 

Nj = 2V^ s Q ^ - — —J cosd. (14) 

Setting (N 0 V sat /s 0 Vj ) = 1 for test purposes one has (13). The relation between 
angle of attack, 6 , and roll angle 0 = cut, is given by 

cos# = acoscot - j3 (15) 

as shown in Figure 3. Also a and (3 are constants determined from satellite ge- 
ometry and the minimum angle of attack. 

It is then possible to integrate G(s) both numerically and analytically, as- 
suming some shift, 0 , in the roll axis: 

f±075 + 5 f±075 + 5 

1+ = I G(s) d0 = 2\/7r s 0 I (acos0-)3)d0 (16) 

J±0 35 + 5 J±0 35 + 5 

where ±0g is determined by inverting (15). Then exploiting G(s)'s symmetry 
about 0-0 and using trig identities one has: 

I + + I_ a [cos0 75 - cos03 5 ] sin 5 

I + - I_ a [sin0 75 - sin0 35 ] cos 6 - ^ (0 75 - 0 35 ) 


as 


1 3 (075 “ 035 ) 
a (sin0 75 - sin0 35 ) 


« 1 


(17-a) 


for SMC-NACE geometry in which the orifice normal is normal to the satellite 
spin axis to within 0.5°, one has : 


\ 

and then: 


5 = 


arctan 


sin0 75 - sin0 35 
COS0 7 5 — COS 035 



(18) 


= (I + - !_)/ {A\fn a (sin0 75 - sin0 35 ) cossj 


( 19 ) 
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As mentioned above the s 0 determination is equivalent, within some scale factor 
(in this case set to unity), to. determi ning density. Hence, density may be calcu- 
lated by this method. 

The shift angle, 8 , determined is the combined shift resulting from the elec- 
trometer transformation plus any errors in the telemetry data. 

In practice the integrals, I±, are computed using the trapezoidal rule: 


I + 


I 


co At 


-co At 


(G(tj) + G(tj + ! )) 

|_t 35 < tj < t 75 

(G(t i )+G(t i _ 1 )) 

-t 35 >ti>-t 7 5 


( 20 ) 


( 21 ) 


where: tj = time measured with respect to the origin (0 = cot = 0) given by the 
telemetry data 


A t = data spacing (in seconds) 
co = 2 7r/T. 


Thus, it can be seen that the bulk of the calculation consists of merely summing 
the data points in the regions of time used. 

In general, one may expect that the initial guess of the origin may be dis- 
placed far enough that the assumed conditions on s and 8 are violated in the inte- 
grations. The calculated shift and density will then be only approximately cor- 
rect. To compensate for this, the technique may be applied iteratively using the 
successive 5’s to adjust the integration limits. 

In testing this procedure a G(s) curve was generated with s = 10.052 and 
At = 0.001 sec, characteristic of the experiment data. This curve was then 
shifted 60° in roll angle and the minimum angle was set to zero to maximize the 
errors induced by integrating over this region. The results are shown in. the 
table in Figure 7. The technique has converged in two iterations after the initial 
integration to accuracies much better than the experimental accuracies, (i.e. a 
residual phase error of ~ 0.01° and a 0.05% error in s 0 . It can be expected that 
in general only 1 iteration will be required as phase shifts greater than 30° should 
be rare . 
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Possible sources of systematic error in this method arise from the term 
X = (3 (4> 1S - 0 35 ) in equation (17-a), i.e. the deviation of the angle between the 
orifice normal and the satellite spin axis from 90° . There is also an error in- 
duced by the variation in ambient density over a spin period. 

In considering the orifice normal-satellite spin axis deviation, the expected 
deviation from the nominal value 0. 5°, was used. This resulted in a zero re- 
sidual phase shift (curve symmetry is preserved) and a density error of 0.37% 
(under estimate). It should be noted that if this deviation of density was large 
(e.g. a non-equatorial orifice mounting) and was known from attitude data, equa- 
tion (19) can be used to correctly find the source density by including the term 
initially dropped from the denominator in (17). 


The variation in density over a spin period is given by: 


exp 


^-28.0 V sa t >n t 


N H 


( 22 ) 


where: t is the time in seconds measured from the spin period's center, H, N 
are the mean density scale height and the mean molecular weight, 28 is the mo 
lecular weight of N 2 , V sa t » n is the vertical satellite speed. Using the nominal 
SM-C orbit the exponential effect was maximized by using: 


V sa t,n = 0.3 km/sec 
H = 37.86 km 

N = 11.79 

Where H and N are taken from the 1971 Jacchia model thermosphere for a 400 km 
altitude and an exospheric temperature of 600° K. Using the nominal satellite 
period one has : 

exp (- 0.0188 t), - 4 < t < 4 sec. 

This effect was found to result in a residual error in phase angle of 0.85 degrees 
with s 0 accurate to 0.05%. The effect of this altitudinal density variation is to 
give a cosine fit of the same amplitude but phase shifted, with respect to a con- 
stant density fit. 

In conclusion, this method allows the rapid determination of both roll phase 
angle and source density to a high degree of accuracy. The density determina- 
tion is performed using a sample of about 50% of the data points in the angle of 
attack range -90 to 90 degrees and also implicitly averages the noise over very 
long intervals. As it provides a. density characteristic of a spin period it s only 
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deviation results from the ambient density variation over a cycle. Since the 
exponent is small the variation is approximately linear: 

exp (-0.0188 t) ~ 1.0 - 0.0188 t 

the assignment of the density to the midpoint of the spin period thus gives the true 
density at the midpoint with errors of order (0.0188t) 2 . In Figure 8, G(s) and a 
cosine density fit are shown. 

One additional error source for density determinations is that due to the un- 
certainty in minimum angle of attack, 6 min, of approximately ±2 or 3 degrees. 
This error enters through the cos(0min) dependence of a in equation (19). For 
example, for0 min = 20° a deviation of ±2% may be expected. This is apparently 
the largest error contribution of the three that were investigated and thus sets 
an upper limit on the accuracy that may be achieved. 


V. NOISE SMOOTHING: INTERVAL AVERAGING 

A second method of spin modulated density determination considered was 
interval averaging. In this technique a simple average: 

GO) = jjj 2 Gj(t) (23) 

is taken over some odd number of points with the average being assigned to the 
interval’s midpoint, t . This process is also known as taking a "running mean" 
of a curve. If the roll angle phase is determined as in Section IV, then it can be 
shown that this method may be used to average out noise which yields source den- 
sities upon division by the spin modulation factor. 

It was desired to find the interval size that would both smooth the data best 
while simultaneously least distorting it. This goal resulted in an investigation 
of the following considerations for the choice of an interval size, 5 1: 

(1) The error induced by the averaging processing itself, 

(2) The degree of noise smoothing obtained, and; 

(3) The amount by which the geometry correction factor is "washed out". 

A computer analysis was made of the variation of fractional error, f as a 
function of angle of attack, where : 

f =(G(t) - G(t)) /G(t) (24) 
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The results, as shown in Figure 9 for an interval size of 0. 1 seconds and data 
sampling rate of 0.01 seconds, were that: 


(1) In the linear s region f was constant, for t} Sq (0) ~ 1, quite small and 
negative, 

(2) For t? So (0) > 1 f showed a characteristic oscillatory structure signifi- 
cantly greater than in (1) but still small, 

(3) In the exponential s region the error is exponentially increasing and 
positive. 

Thus, f essentially measures the deviation of G(s) from linearity in s. In the 
exponential s region there exists a systematic error in that the linear average 
of an exponential function will always be greater than the function evaluated at 
the interval’s midpoint whenever the function is monotonically increasing or de- 
creasing over the interval. 

As it is only the linear s region that is of interest in density determination, 
and as it is only here that the averaging process does not contain the aforemen- 
tioned systematic error, only for this region will a quantative relation between 
f and 5t be derived. 

If one now considers a Taylor series expansion of f about the interval's mid- 
point for 


G(t) = 2>/7r s 0 a cosco (t + tj) 


(25) 


by symmetry, odd powers of cot will cancel when summed over the interval and; 


_ -co 2 y-v tj 2 

f Ai " 2N ^ — 5t/2 < tj < St/2 


(26) 


where : 

0 0 - satellite spin frequency (rad/sec) 

N = number of points in averaging interval 
5t = averaging interval size 

tj = time of the ith data point measured from the interval's center (sec) 
At = data point spacing in seconds. 
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Note that this result is independent of s. Figure 10 compares the values of fwith 
those calculated by computer. The result of neglecting the remaining terms of 
order (cotj) 4 is an error in the predicted f of ~10 -4 . 

As can be seen in the table of Figure 10, the formula is of high precision. 
This formula can therefore be used to estimate the effects of data compression, 
i.e. not averaging every point but selecting out some subset. For example, sup- 
pose 1 out of every 10 data points was chosen for averaging then the increase in 
error is given by the ratio f Q i /f O oi * as a f unc ^ on averaging interval 
size which can be calculated. The result of this calculation is shown in Figure 
11. As the error increase is small, preliminary selection may be desired to 
minimize computation time. 

Consider the final case of tj Sq ( 6 ) > 1 for small angles of attack. Retaining 
the linear and cubic terms in the expansion, one has: 


where : 


f At =j - — — V* 7?s (coCt + tj)) [1.0 -(tancot) cotj 

/ Nt? Sq (cot) ° 

-(cotj) 2 /2 + (tancot) (cotj) 3 /6] i - 1. 



< 


5t 

2 


(27) 


Note that f now depends on both s Q and roll angle. In Figure 12, f 01 and f 0Q1 
are plotted vs. angle of attack for 5 t = 0.1 sec. For f 01 the result agreed 
with the computer calculation in Figure 9 to within 1%. It is also notable that 
the increase of f with increasing At (data selection) is small. 


Figure 12 also demonstrates the "washing out" effects in the smoothing down 
of the G(s) peak (See Figure 4) and the enhancement of the lower sides of the cor- 
rection as it is smeared out more uniformly by the averaging process. 


In summary, from a consideration of systematic averaging errors, the most 
severe is that produced by the "washing out" of the geometry correction if this 
is to be studied. 


In the actual noise smoothing itself, it is reasonable to ask what is the best 
interval size in terms of minimizing the noise contribution to the source density. 
As shown in Figure 2, it would appear that by inspection most of the noise am- 
plitude is contained in the higher harmonics. Given this fact, consider the am- 
plitude of the Nth noise harmonic, H N of a Fourier series representing the post- 
electrometer noise. 
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Where: G (s) G (s) + H (s) 

with noise ideal noise’ 

Then H N = a N • cos (cot+ </> N ), and for a fixed averaging interval size, 6t, the 
greatest error magnitude contribution to the mean is: 

a " ; h,lf cycle COs( “' + *N> 

Tt < 28 > 

Where: a N is the constant coefficient of the Nth harmonic, and, 

0 N is the phase angle of the Nth harmonic. 

Thus, the contribution is minimized by making 5t as large as possible. 

In conclusion, for interval averaging as a method of minimizing noise con- 
tributions to source densities, there are two separate regions of interest: 

(1) For rjs 0 (0) ~ 1 (|0| > 35°) the only restraint is on the magnitude of f. 
This should allow interval sizes as large as 0.5 sec with systematic 
errors of <1%. 

(2) For rj So (0) >1, (|0| < 35°) and electrometer bandwidths >8 Hz, the 
most important limitation on interval size is the geometry correction 
being "washed out". This appears to set an upper limit of approximately 
0.1 sec. 


CONCLUSION 

The electrometer transformation has been studied and found to be well ap- 
proximated by linear transformation in which only the roll angle is shifted for 
angles of attack 6 , such that 75 > |0| > 35 degrees. This also applies for \ 6 1 
< 35 degrees when the electrometer bandwidth is > 8 Hz. For small bandwidths 
(e.g. 1.3 Hz) the transformation is non-linear for |0| < 35° and must be treated 
separately. 

Exploiting the fact that the essential cosine curve representing G(s) for 
35 < |0 1 < 75 degrees when passed through the SM-C NACE electrometers is 
preserved, an iterative integration scheme was developed to fit G(s) with a co- 
sine curve. This method rapidly and accurately determines the roll angle phase 
and a representative source density for the spin period. This scheme also im- 
plicitly greatly dampens noise contributions to density for reasons given in sec- 
tion V, as the averaging intervals (i.e. the intervals of integration) are on the 
order of 0. 8 sec. 


12 



Finally, interval averaging was studied as a possible means of examining 
the geometry correction effect for large bandwidth electrometers (> 8 Hz). The 
restraint here was primarily the washing out of detail that limited the interval 
size to ~0.1 sec for the SM-C nominal spin period of 8 sec. 
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Figure 2. San Marco-C NACE Data Channel 7, Molecular Nitrogen Showing Noise Structure 
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Figure 3. Satellite Geometry 
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Figure 6. Effect of Shifting G (s), Post Electrometer, by Phase Shift in Fundamental Harmonic 
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0 

-34.9801° 

34.9801° 

25.0199° 

12.3838 

.231974 

1 

-23.9264 

58.9065 

1.0935 

10.4299 

.037594 
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- 1 .0821 

59.9886 

.0114 

10.0576 

.000557 
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- 1 .26 x I0' 5 

59.9886 

.0114 

10.0573 
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- 1 .26 x I0‘ 5 

59.9886 

.0114 

10.0573 

.000527 


Figure 7. Sample of Integration Results using Theoretical G (s) s 0 = 10.0520, A=0, 0 min =0, 5 $hift = 60° 
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(G (POST)— G (PRE) )/G (PRE) 


COMPARISON OF INPUT G(S) AND So» 2. SQRT (PI) » COSINE (CUi T) 
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Figure 8. Density Determination and Phase Shift Calculation 
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Figure 9. Fractional Error vs. Angle of Attack 
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FRACTIONAL ERROR DIFFERENCE 
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200 

-1.131 x IO" 3 
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.0015 


Figure 10. Comparison of Computed and Predicted values of f for Region I 
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Figure 11. The Effect of Sampling on Fractional Error 
for Interval Averaging 
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Figure 12. The Effect of Sampling on 
Fractional Error for Interval Averaging 
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